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Abstract. The effects of adsorbates on nucleation and growth of two- 
dimensional islands is investigated by kinetic Monte Carlo simulations and 
rate equation theory. The variation of island morphology with adsorbate 
parameters is discussed and the temperature-dependence of island density 
in the case of immobile adsorbates is studied in detail. A set of rate equa- 
tions for the description of nucleation in the presence of predeposited mobile 
and immobile adsorbates is developed. 



1. Introduction 

While a satisfactory understanding of nucleation and growth in homoepi- 
taxial systems is beginning to emerge [1], the dramatic effects of tiny traces 
of adsorbates continue to surprise researchers. Among a wealth of recent 
examples, we may mention the influence of CO on island shapes and inter- 
layer transport on Pt(lll) [2], the modification of attachment kinetics by 
H preadsorbed on Si(OOl) [3] and the H-induced enhancement of Pt self- 
diffusion on Pt(llO) [4]. In this situation it seems useful to explore different 
generic scenarios for adsorbate effects on submonolayer growth within a rea- 
sonably simple, yet flexible, theoretical model. The present paper reports 
on an ongoing effort aimed in this direction. 
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In Section 2 we briefly review our earlier findings and present new re- 
sults concerning the dependence of the island density and morphology on 
the temperature and the characteristics of the adsorbate (mobility and the 
strength of interaction with adatoms). In the Sec. 3 we sketch a rate equa- 
tion theory which provides a unified description of two-dimensional nucle- 
ation in the presence of predeposited mobile and immobile adsorbates, and 
derive analytic solutions for some special cases. 

2. Kinetic Monte Carlo Simulations 

2.1. MODEL 

We employ a recently introduced solid-on-solid growth model with two sur- 
face species A and B, representing the adatoms of the growing material and 
the adsorbates, respectively [5, 6]. The simulation starts on a flat substrate 
composed only of A-atoms. The basic microscopic processes are deposition 
and migration; desorption is not allowed. Two deposition modes can be 
studied: (i) simultaneous deposition (codeposition) of both species, with 
generally different fluxes, Fa and Fb, and (ii) predeposition of a certain 
adsorbate coverage 6b prior to growth. 

Adatoms and adsorbates migrate with a nearest-neighbor hopping rate 
Rd = koe~ ED / kBT , where ko = 10 13 Hz and the energy barrier for a particle 
of type X = A or B is given by 

E%= E (^E^ + n^E^). (1) 

Y=A,B 

Here E^£ is the hopping barrier for a free X atom on a substrate atom Y, 
tiq is equal to one if the substrate atom is of type Y and zero otherwise, n* Y 
is the number of nearest-neighbor X-Y pairs, and E^ Y is the corresponding 
contribution to the barrier (symmetric in X and 1"). There are no lateral 
interactions between adsorbate atoms {E EB = 0). The diffusion rate of a 
free adatom is D = koe~ E ^b/ k B T '. 

In previous work [5], we have shown that in addition a process of place 
exchange between an adsorbate and an adatom with at most one lateral 
bond is necessary to achieve decorated island edges. This occurs at a rate 
koe~ Ee */ kBT . When the barrier E ex is sufficiently low, then adsorbates are 
floating on the edges of growing islands. 

Adsorbates can be mobile or immobile and strongly or weakly interact- 
ing with the adatoms. Their mobility is controlled by the strength of inter- 
action of adsorbates with the substrate which is described by the energy 
barrier E^. The interaction with adatoms is determined by the energy 
barrier ■ 
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2.2. SUMMARY OF PREVIOUS RESULTS 

We have investigated previously the case of floating mobile adsorbates in 
the situation when the concentration of adsorbates is comparable with that 
of the growing material [6] . We studied the dependence of the island density 
iV on the flux Fa and the coverage 9 a for both codeposition and predepo- 
sition. We found that the adsorbates strongly increase the island density 
without appreciably changing its power-law dependence on flux, the in- 
crease being stronger for larger E^ B . The increase was only slightly higher 
for predeposition than for codeposition. 

We also observed a stronger coverage dependence of the island density 
in comparison with homoepitaxy. This was interpreted as a delay of the sat- 
uration regime, where island density becomes independent of coverage. The 
coverage dependence is again more pronounced for strongly interacting ad- 
sorbates, but it is weaker in the case of predeposition than for codeposition. 
A further noteworthy feature is the much higher density of free adatoms 
than in homoepitaxy, which shows an intriguing, oscillatory coverage de- 
pendence in the case of codeposition [6]. 

Many of these simulation results could be qualitatively explained in 
terms of a simple rate equation theory, which rests on the assumption that 
the adsorbates affect the growth process only by slowing down the diffusion 
of adatoms. If, in addition, the adsorbates are treated as immobile traps 
(which may be justified at least for sufficiently high adsorbate coverage 9b, 
see Section 3.4), the adatom migration is described by an effective diffusion 
coefficient 

° {9b) = i-e B + e B e^ B i^T- & 

Recently, we have considered a situation when only a low concentration 
of predeposited adsorbates is present on the surface [7]. We investigated 
the dependence of island density on flux and on the concentration of ad- 
sorbates 9b- Even concentrations as low as 9b = 0.002 ML can lead to 
a severalfold increase of the island density. However, the increase signifi- 
cantly depends on the mobility of the adsorbates. In the case of essentially 
immobile adsorbates = 5 eV), a new feature in the flux dependence 

of the island density appears. Instead of a single power law relationship, 
there is a plateau where N ~ 9b, reflecting the dominance of heteroge- 
neous nucleation in a flux interval F\ > F > F<i- Scaling arguments yield 
the estimates F\ k, D9 b , F2 ~ 9 B De~ En / fcflT for the plateau boundaries 
in terms of impurity coverage and strength. For sufficiently small 9b these 
were confirmed by scaling plots of N/9b vs. F/F\ and F/F2 [7]. 

In this work we further extend the above results. In the following sub- 
section, we discuss effects of adsorbates on the island morphology, and in 
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the subsequent subsection we present results on the temperature depen- 
dence of the island density. The energy barriers = 0.8 eV, E^ = 0.1 
eV, E BB = o.l eV, and E AA = 0.3 eV were fixed in the following. Other 
energy barriers as well as the ratio D/F were varied. 

2.3. ISLAND MORPHOLOGY 

The island shape in homoepitaxy depends on the rates of the various kinetic 
processes involved. The size of islands varies with the ratio D/F. To be 
specific, we set in this subsection the temperature to T = 500 K. Then, 
the island shapes in homoepitaxy with the parameters mentioned above in 
Section 2.2 are regular (approximately rectangular, close to square shape, 
nondendritic) for fluxes in the interval F = 0.00025 — 0.64 ML/s. The 
shapes become irregular for higher fluxes. In order to evaluate the effect 
of adsorbates, we performed several sets of simulation with fixed flux F = 
0.001 ML/s and varying adsorbate properties. 

Let us consider first the case of floating adsorbates. It is useful to dis- 
tinguish two situations: (i) well decorated island and (ii) poorly decorated 
islands. In the former case, the majority of edge sites is decorated by an 
adsorbate. We found that the shape is very regular, only the islands are 
smaller; in fact it can be shown that the decoration enhances edge diffu- 
sion and thus makes the island edges smoother [6]. The shape is not very 
sensitive to the mobility of adsorbates. It only has to be sufficiently high 
so that adsorbates can reach the island edges. The exchange of adatoms 
approaching the step edge and the adsorbates decorating the step edge is 
the dominant mechanism. 

If the islands are poorly decorated (the number of adsorbates attached 
to the island edges is considerably less than the number of perimeter sites) , 
then, the shape changes to irregular and regularity is restored again only 
for very small 0b- Rather than promoting edge diffusion, in this regime 
the adsorbates bound on the island edges act as sinks for adatoms provided 
that the adatom-adsorbate interaction (given by E AB ) is sufficiently strong. 
Another factor affecting the island shape in this situation is the change of 
the mobility of adsorbates (given by Ef^). Examples of the different island 
mophologies for various values of parameters E AB and E B ^ are shown in 
Fig. 1. We can see that the increase of both energy barriers leads to more 
irregular island shapes and also to an increase of island density. 
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Figure 1. Examples of configurations obtained with flux Fa = 0.001 ML/s at a coverage 
9a =0.1 ML after predeposition 9 B = 0.005 ML of adsorbate coverage, E AA = 0.3 eV, 
iJ ox = 1 eV and different energy barriers E AB — 0.2 eV (left column), E AB = 0.4 eV 
(middle column), E AB = 0.6 eV (right column), E B U A : from top to bottom Ef u A = 0.7 eV, 
E B U A = 1.0 eV, E B U A = 1.2 eV, E B U A = 1.5 eV. We show only 50 x 50 sections of larger 
simulation boxes. 
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The irregularity is not so apparent for weak adsorbate-adatom interac- 
tion (left column) or high mobility of adsorbates (top row). The immobile 
strongly interacting adsorbate catches more diffusing adatoms and at the 
same time it is floating on the growing island edge. Hence, the island grows 
faster in the neighborhood of such adsorbates. This kind of adsorbates 
also block the motion of adatoms along the edge and thus suppresses the 
smoothening of the island edges. 

Let us discuss briefly the island morphology for non-floating adsorbates. 
The mechanism of floating becomes less relevant for adsorbates with low 
mobility and does not act at all in the limit of immobile adsorbates. Non- 
floating adsorbates are incorporated inside a growing island and beyond a 
certain island size have no effect on growth. In this context, we want to note 
that surface defects can be treated in our approach as immobile adsorbates. 

2.4. TEMPERATURE DEPENDENT ISLAND DENSITY 

While in our simulations we have so far focused on the flux dependence 
of the island density [5-7] in experiments it is most often measured as a 
function of temperature at fixed flux. Therefore, we complement our previ- 
ous results by the calculation of the temperature dependence. We restrict 
ourselves here to the case of predeposition of immobile (non-floating) ad- 
sorbates. In Fig. 2 we show data for several adsorbate concentrations. We 
find, similar to the flux dependence [7], a plateau interval T\ < T < T2, 
where N is almost independent of T. For low temperatures, T < T±, ho- 
mogeneous nucleation is more important than heterogeneous nucleation by 
adsorbates, because the adatoms are not sufficiently mobile to reach the ad- 
sorbates. In the plateau region nucleation is predominantly heterogeneous. 
For high temperatures, T > T2, adatoms start to detach from adsorbates 
and homogeneous nucleation is again important. This kind of behavior is 
observed experimentally in defect nucleation on oxide and halide surfaces 
[8] and in growth on surfaces with strain-relief patterns, which act as an 
ordered array of defects [9] . 

For a quantitative estimate of the transition temperatures T\ and T2, 
we recall the three relevant time scales in the problem [7] : (i) The diffusion 
time td ~ 1/D9b, which is the time required for an atom to explore the 
"capture zone" of area 1/0b associated with a single impurity, and hence to 
get trapped at the impurity, (ii) The deposition time tf ~ &b/F, which is 
the time between subsequent arrivals of adatoms within the capture zone, 
(iii) The trapping time r ~ (l/D)e En / fcsT at the impurity. The plateau 
regime is characterized by r > tf > td, which yields the expressions 
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Using = 0.8 eV, F = 0.001 ML/s and 9 B = 0.001 ML, we get T x = 
403 K, while T 2 = 465 K for E^ B = 0.4 eV and T 2 = 542 K for B = 0.6 
eV, in good agreement with the simulations. From the expression for T 2 , 
we directly see how the width of the plateau grows with increasing E^ B ■ 




1.86 1.96 2.06 2.16 2.26 2.36 2.46 

1000/T 

Figure 2. Averaged island density as a function of the inverse temperature for different 
concentration 9 b of predeposited impurities and different energy of interaction between 
adsorbates and adatoms: E^ B = 0.4 eV - filled symbols, E£ b = 0.6 eV - open symbols. 
The adatom interaction energy E^ A = 0.3 eV and the energy barriers E B U A = E cx = 5 eV 
are fixed. Deposition flux is Fa = 0.001 ML/s, the coverage of adatoms is 6a = 0.1 ML. 
The behavior in the absence of impurities (homoepitaxy, 9b = ML) is shown for 
comparison. 



3. Rate Equations 

In this section we develop a set of rate equations for the description of two- 
dimensional nucleation in the presence of mobile and immobile adsorbates. 
Throughout, the impurities will be assumed to have been predeposited at 
a coverage 9b, though codeposition is easily incorporated along the lines 
of our earlier work [6]. To keep matters transparent, we also restrict our- 
selves to the case of irreversible aggregation (critical nucleus size i* = 1). 
Detachment of impurities from adatom islands will, however, be included, 
since this process plays a crucial role at low deposition fluxes. 

In addition to the adatom density n and the island density N, a minimal 
description of nucleation in the presence of adsorbates turns out to require 
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the monitoring of two additional quantities [8, 10]: The density of free 
impurities no, and the density n c of critical clusters bound at impurities; in 
the present case i* = 1 the latter is simply the density of bound adatom- 
impurity pairs. Taking into account all relevant processes, one arrives at 



the following set of evolution equations: 

dn/dt = F - Dn(2n + N) - Dn(n + n c ) + n c /r, (4) 

dN/dt = Dn(n + n c ), (5) 

dn /dt = -Dnn - Dn N + n c /r + N/t', (6) 

dn c /dt = Dn{riQ — n c ) — n c /r. (7) 



Here 1/r ~ kQe~^ E <™b +E ™ B )/ k BT d eno t e s the rate of detachment of an 
adatom from an impurity, D = koe~ E eub/ kBT is the impurity diffusion con- 
stant, and 1/r' « fc e~ < - £ ' sub+£ ' n )/ fcsT is the detachment rate of impurities 
from islands. In the expression for r, we have assumed that the impurities 
are (much) less mobile than the adatoms, so that the relative mobility is 
dominated by the adatom mobility. 

Most of the impurity terms in Eqs.(4-7) should be self-explanatory. We 
merely point out the term n c /r in (4,6,7) which describes the dissociation 
of adatom-impurity pairs, and the term N/t' in (7) which describes detach- 
ment of impurities from the islands. Strictly speaking, all reaction terms 
should be adorned with dimensionless capture coefficients. However, since 
our objective here is a qualitative, rather than a quantitative description, 
all these coefficients have been set to unity. 

3.1. DIMENSIONLESS FORMULATION 

Following Tang [11] we now rescale the adatom and island densities by 
y/F/D and time by 1/ V DF. In addition, the impurity densities no and n c 
will be rescaled by the initial impurity coverage Ob- Marking dimensionless 
quantities with a hat, we obtain the rescaled evolution equations 



dn/dt = 1 — h(2n + TV) — an(fiQ + n c ) + bh c , (8) 

dN/di = n(n + an c ), (9) 

dho/di = -nh - cn N + (b/a)h c + (bc/a 2 )N, (10) 

dh c /dt = h(ho — h c ) — (b/a)h c , (11) 



where a = 9b\/D/F, b = 6b /Ft and c = D/D. The initial conditions are 
n(0) = N(0) = n c (0) = 0, n (0) = 1. 
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The parameters a and b may be written in terms of the characteristic 
fluxes, F\ and F2, which limit the plateau regime in the case of immobile 
impurities, as a = y/Fi/F and b = F2/F. For impurities to be at all 
relevant, it is necessary that F1/F2 = D9bt = a?/b » 1. By retaining 
only the most important terms, we can analytically extract the behavior 
predicted by Eqs.(8-ll) in simple cases. So far this has been achieved only 
for immobile impurities (c = 0). 

3.2. IMMOBILE IMPURITIES WITHOUT DETACHMENT 

Here we consider the plateau regime, where a» 1 and 6^1. Numerical 
integration of (8-11) with b = c = and a » 1 shows that the adatom 
density is time independent, taking the value n = F/D9b to high accuracy. 
For early times, this reflects the balance between deposition and capture 
at impurities, while at late times the capture at islands dominates; since 
N = 9b at late times, the resulting adatom density is the same in both 
regimes. Setting n = F/D9b, the remaining three rate equations become 
linear and are readily integrated, with the result 



Both impurity species decay exponentially at late times. 
3.3. IMMOBILE IMPURITIES WITH DETACHMENT 

The regime F <C F2 is detachment-dominated, in the sense that an adatom 
typically visits many impurities before being incorporated in an island. 
The primary effect of the impurities is then to slow down the diffusion by 
temporarily trapping the adatoms [6]. To see how this emerges from the 
rate equations, we infer from numerical integration that, in the relevant 
late time regime, an equilibrium is established between attachment and 
detachment of adatoms at free impurities. This implies that Dnno « n c /r. 
Equation (6) then yields dno/dt ~ 0, so that the impurity concentration 
essentially retains its initial value no = 9b, and n c ~ (D9bt)ti. In this way, 
the impurity rate equations (6,7) have effectively been eliminated, and we 
are left with a pair of modified equations for n and N. Taking into account 
that D9bt = F1/F2 3> 1, and therefore n c » n, they read 



N = O B (l-(l + O/O B )e- ' 0B ) 

n = 9 B e- e/(>B , 
n c = 9e- e ^. 



(12) 



(13) 
(14) 



dn/dt = F — Dn[(D9 B T)n + N] 
dN/dt = (D9 B T)Dn 2 . 



(15) 



(16) 
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Standard analysis [11] shows that the asymptotic island density is of the 
order N ~ (FtOb) 1 ^ 3 , which can be brought into the familiar form N ~ 
(F/D) 1 / 3 by identifying the effective diffusion constant as D = 1/(t#b). 
This is precisely the strong impurity limit of (2). 




Figure 3. Island density as a function of deposition flux, obtained by numerical solution 
of the rate equations (4-7) for different choices of the impurity diffusion barrier. The 
adatom coverage is 9 = 0.1 ML and the impurity coverage 6b = 0.005 ML. 



3.4. THE EFFECT OF IMPURITY MOBILITY 

The second term on the right-hand side of Equation (6) describes the de- 
pletion of free impurities due to capture at stable islands. This is partly 
compensated by the detachment of impurities from islands described by 
the last term, but nevertheless, the net effect is expected to be a decrease 
of the island density, because the depletion reduces the ability of mobile 
impurities to act as nucleation centers. 

In the absence of an analytic study of the full set of rate equations, in 
Fig. 3 we provide some sample results obtained by numerical integration. 
The adatom diffusion constant D and the detachment time r were chosen in 
accordance with the parameters employed in the KMC simulations. Data for 
three different choices of the impurity diffusion barrier E^ are shown, and 
the behavior in the absence of impurities is included for reference. The case 
of low impurity mobility (E^£ = 1.2 eV, corresponding to D/D rs 10 -4 ) 
is particularly noteworthy, as it indicates an intermediate scaling regime 
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where the island density scaling exponent, defined by N ~ F x , is larger 
than the homoepitaxial value \ = 1/3 (a fit yields x ~ 0.42). Further 
exploration of this phenomenon will be presented elsewhere. 
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